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Abstract 

We demonstrate thermo-mechanical relaxation oscillations in a strongly driven quartz crystal. 
Dynamic bifurcation leads to two stable oscillation states with a distinct electrical impedance. 
Slow Joule-heating, which shifts the susceptibility of the crystal, provides a feedback that leads to 
thermally-induced oscillations, in which the amplitude of the crystal is modulated by a relaxation 
cycle. The frequency of the relaxation cycle is roughly a million times lower than the resonance 
frequency of the crystal, and it can be adjusted by the detuning from the critical point for dynamic 
bifurcation. The experimental observations are reproduced by a simple model that takes into 
account the slow dynamics of the system. 
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Harmonic oscillators, such as quartz crystals, are used for timing and sensing purposes, 
and constitute an indispensable part of modern electronic devices. Besides harmonic oscil¬ 
lators, other types of oscillators such as the relaxation oscillator occur frequently |T]. Where 
harmonic oscillators conserve energy — they exchange kinetic with potential energy, as in a 
mass-spring system- a relaxation oscillator repetitively dissipates energy. It’s dynamic be¬ 
haviour can be described by coupled hrst-order differential equations, with non-oscillatory 
steady state limits. The period of a relaxation oscillator is determined by energy decay 
rates, which result, for example, from viscous friction or resistive dissipation. 

Canonical examples of relaxation oscillators include the electronic flip-flop [2], the 
mammal heart muscle [3], and the Pearson-Anson neon-lamp oscillator |1]. Mechanical 
relaxation oscillators have been observed recently in a micro-electromechanical impact de¬ 
vice [5], a held-emitting carbon nanotube |6], and in a system of coalescing nanofluidic 
droplets [7]. Here, we report a dynamic relaxation oscillator, with oscillatory steady states, 
in which the oscillation amplitude of a nonlinear resonator is modulated by a relaxation 
cycle. The period of the relaxation cycle is determined by the thermal time constant and 
the mechanical ring-down of the device. 

The dynamic relaxation oscillator is implemented using a commercial AT-cut quartz 
crystal with a specihed resonance frequency of 4.608 MHz at the fundamental thickness 
shear-mode. The crystal is driven by an rf-voltage generated by a lock-in ampliher, and 
its motion is transduced using a current probe |H] and detected by the lock-in, as is shown 
in Fig. 1(a). The experiments are performed at room temperature and in atmospheric 
pressure |9]. Figure 1(b) shows the response of the weakly driven quartz crystal, which 
corresponds to a harmonic resonator, with a resonance frequency of fo = 4.607 MHz and a 
mechanical ring-down time of Tm = QQ^vrfo ~ 50 ms. 

When the driving voltage is increased, the response starts to deviate from a harmonic 
oscillator, as the resonance peak becomes non-symmetric. The resonance frequency shifts to 
a higher value: a signature of nonlinear behaviour with a positive higher order spring con¬ 
stant. The anharmonic behaviour of a quartz crystal has been studied previously mm, 
and in a recent quartz crystal microbalance (QCM) experiment, a nonlinear response was 
used to enhance the responsivity to an added mass [12]. 

When the crystal is driven beyond the critical point, at Vref = 0.02 V and 
fref = 4 .6071 MHz, the nonlinearity gives rise to a bifurcation and hysteresis occurs: 
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two stable states co-exist, and the crystal oscillates either at a high or at a low amplitude. 
Figure 1(c) shows a hysteretic frequency response of the quartz crystal. When the frequency 
is swept from a low to a high value, the crystal susceptibility follows the upper branch, and 
it oscillates at a high amplitude. For a reverse sweep, the low amplitude state is stable. 
Figure 1(d) shows the driving conditions that give rise to a bistable response. The hgure 
is constructed by superimposing frequency response measurements taken in forward and in 
reverse direction (horizontal, fast axis), at a varying drive amplitude (vertical, slow axis). 
The hysteretic transitions between the two oscillating states, which occur at the bifurcation 
frequencies f I and f form the switching element that is required for the relaxation 
oscillator. 

The electrical impedance of a quartz crystal depends on the amplitude of the 
oscillation, and it is bistable when the amplitude is bistable, as in the regime visualized 
in Fig.l(d). We measured that Ziow~ 80011 in the low state and Zhigh~20r2 in the 
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FIG. 1: (a) Schematic of the measurement circuit; inset: photograph of the quartz crystal, (b) 
Driven response (amplitude shown only) at the fundamental shear-mode driven at Vref = 0.01 V. 
The arrows indicate the sweep direction, (c) Bistable response when driven at 0.3 V. (d) Drive 
parameters that lead to a bistable shear-mode oscillation. The figure is obtained by superimposing 
forward and backward traces, such as the one shown in (c) (horizontal, fast axis), while varying 
the drive power (vertical, slow axis). Color scale: susceptibility of the quartz crystal. 
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FIG. 2: (a) When the crystal is driven close to the bifurcation point inside the hysteresis regime, 
it decays to the low-amplitude state, (b) When driven close to the bifurcation point outside the 
hysteresis regime, relaxation oscillations occur. 

high state. Given the constant driving voltage, the dissipated power then depends on 
the oscillation state. Joule-heating of the crystal gives rise to a signihcant change in its 
mechanical properties. In particular, the negative temperature dependence of the Youngs 
modulus of silicon oxide causes the mechanical spring constant, and thus the (nonlinear) 
resonance frequency, to decrease with temperature. As a result, the susceptibility of the 
crystal at the driving tone is increased, and this thermo-mechanical coupling presents the 
feedback that gives rise to relaxation oscillations. 

Figure 2 shows time traces of the crystal amplitude when driven at Vref = 0.5 V at 
two different hxed frequencies. When driven close to f t at the high-amplitude branch, as 
shown in panel(a), the crystal temperature increases, causing the resonance frequency to 
decrease and the susceptibility to increase. The response then follows the upper branch 
as indicated by the arrows, until f ),< hef, where the impedance becomes low. The crystal 
then cools and follows the lower branch until it reaches a steady low-amplitude state. 
Figure 2(b), right hand panel, shows a measured time trace of the crystal amplitude, with 
a single transition to the low state at t = 2s. A similar process occurs when the crystal 
is driven close to ft but outside the hysteretic regime, as in panel (b), but in this case 
the temperature further decreases, until f t> hef forces an upwards transition. Here, the 
process repeats, and the crystal enters a relaxation oscillation. Figure 2(b), right panel, 
shows a time trace of the oscillation, with a period of approximately 1 s. 
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To investigate the relaxation oscillator in more detail, we measured time series of 
the crystal susceptibility close to the bifurcation point, f for a range of driving voltages 
and frequencies. Figure 3(a), shows the frequency of the relaxation oscillation on the color 
scale, for a range of drive parameters. The bifurcation points, as measured in Fig.1(d), 
are indicated by the solid blue lines. In the dark blue area the frequency is zero and the 
relaxation oscillations are absent. Relaxation oscillations occur in a distinct regime close 
to f t, at a frequency that can be adjusted over f^o = 0..2.5Hz by adjusting the drive 
parameters. 
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FIG. 3: Frequency of the relaxation oscillations (color scale) for a range of driving parameters, as 
obtained from experiment (a) and model (b). For each pixel, one time trace is measured from which 
the frequency is determined. In the model, the drive strength, aj, and the frequency detuning, Af, 
are normalized to the critical drive strength and frequency, acr and hr, respectively. 

To corroborate the experimental results, a simple model is developed. While the 
power and frequency of the bifurcation points in a nonlinear resonator can be modelled 
following the analysis by Lifshitz and Cross ra. here we simplify the problem by taking 
only the slowest dynamics into account. We assume a constant susceptibility (dissipation) 
on the high and on the low branch of the hysteresis regime, and a linear dependence of the 
bifurcation frequencies on the drive strength, ft, fj, cc Vref. We take a linear dependence 
of the bifurcation frequency on the crystal temperature, ft, ft oc T [H], and model the 
behaviour using a hrst order heat equation. 
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With the dissipation Plow in the lower branch and Phigh in the upper branch [15], the 
temperature of the crystal is described as [121 



( 1 ) 


where t denotes time, r is the thermal time constant of the crystal, and ATj is the initial 
temperature difference with respect to room temperature. When K is the thermal con¬ 
ductance between the crystal and the environment. Thigh = Phigh/K and Tiow = Piow/K 
represent the two thermal equilibria to which the oscillator relaxes. By solving the heat 
equation and inserting f oc T, one obtains the thermal oscillation period: 



The relaxation oscillation period is thus set by the difference between the thermal 
equilibrium frequencies how and fhigh, and the bifurcation points ft and fj,. Reducing the 


denominator in either of the terms in Eq.2 increases the cooling or heating time, and thus 


the oscillation period, as is observed in the experiment of Fig. 3(a). Equation (2) was solved 
for a range of excitation frequency and amplitudes with t = 1, K = 1, Phigh/Plow = 40, and 
the frequency of the relaxation oscillation, f^O; is plotted in Fig. 3(b). As in the experi¬ 
ment, relaxation oscillations occur for hef < f t- Moreover, the frequency of the relaxation 
oscillation increases with the detuning and it exhibits a maximum. The model could be 
rehned to capture the experimentally observed features in more detail. For instance, the 
frequency-dependence of the temperature in the high and in the low amplitude state could 
be taken into account, and the kink in the f I branch could be modeled by accounting for 
higher-order nonlinearities or multiple (internal) resonance modes [T2] . 

The dissipated power can be estimated from the impedances in the high and low states, 
and we calculate that Phigh ~ 2.5 W and Plow ~ 0.1 W, compared to a dissipation of 1 mW 
in the linear regime of Fig. 1(b). It is interesting to obtain an indication of the temperature 
range during the relaxation cycle. To this end, the crystal was removed from its package, 
and its temperature was measured using an infrared detector mi, facing the crystal at a 
distance of 2 mm. Figure 4 shows the temperature of the crystal during the oscillation. 
When driven at hef = 4.60727 MHz and Vref = 0.55 V in an ambient temperature of 20 °C, 
the mean temperature is 35 °C, and the peak-to-peak value is 7°C. The frequency of the 
relaxation oscillation is slightly shifted compared to a packaged crystal due to a different 
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thermal conductance, and the oscillation is less stable due to the exposure to the environ¬ 
ment. 

Relaxation oscillators are very sensitive to their environment, as is indicated by the 
frequency fluctuations that can be observed in the measurements of Fig.2(b) and Fig.4, in 
which no special precautions were taken for stabilization. This sensitivity could be exploited 
in detectors ca- the relaxation rates respond, besides to mass, to changes in the viscos¬ 
ity and the thermal conductivity of the environment, enabling application as a viscosity 
or a pressure sensor. Since the output of the oscillator resembles a frequency- and pulse- 
width modulated binary signal, the interface to a digital circuit could be simplihed. The 
quartz crystal may also be used as an experimental platform for more fundamental studies 
on dynamic relaxation oscillations, such as the complex but slow behaviour that occurs in 
electrochemical systems [18]. Here, the mechanical device presents interesting dynamics on 
an experimentally convenient time scale. Finally, we note that the thermal time constant 
scales linearly with dimension (volume-to-surface ratio), and that for micrometer-sized de¬ 
vices similar processes could occur at frequencies in the kHz range. 

In conclusion, we demonstrate dynamic relaxation oscillations in a strongly driven 



FIG. 4: Temperature of the quartz crystal (red trace, left axis) during the relaxation cycle (blue 
trace, right axis), measured after removing the crystal from the package. The crystal is driven at 
hef = 4.60727 MHz, and the driving amplitude is Vref = 0.55 V. 

quartz crystal. The amplitude of the crystal oscillates at a frequency that is determined by 
the thermal relaxation time of the crystal, six orders of magnitude lower the fundamental 
harmonic oscillation frequency of the crystal. The dynamic behaviour is captured by a model 
that takes into account the slow dynamics. The frequency of the relaxation oscillation is 
very sensitive to the susceptibility of the crystal, which depends on its mechanical properties. 
This feature may be deployed in a low-frequency readout scheme for QCM-based sensors. 
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